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QQ ' Big Bang nucleosynthesis requires a fine balance between equations of state for photons and 

, relativistic fermions. Several corrections to equation of state parameters arise from classical and 

■ quantum physics, which are derived here from a canonical perspective. In particular, loop quantum 

' gravity allows one to compute quantum gravity corrections for Maxwell and Dirac fields. Although 

the classical actions are very different, quantum corrections to the equation of state are remark- 
ed ■ ably similar. To lowest order, these corrections take the form of an overall expansion-dependent 
multiplicative factor in the total density. We use these results, along with the predictions of Big 
Bang nucleosynthesis, to place bounds on these corrections and especially the patch size of discrete 
. quantum gravity states. 

PACS numbers: 04.20.Fy, 04.60.Pp,98.80.Ft 

^' 

6 

^ , Much of cosmology is well-described by a space-time near a spatially isotropic Friedmann-Robertson- Walker models 
22 ' with line elements 



■ 
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I. INTRODUCTION 



ds^ = -dr^ + a(T)2 ( + {M^ + sin^ dd^^) ] (1) 



> , 

: where k = or ±1, sourced by perfect fluids with equations of state P = wp. Such an equation of state relates 
' the matter pressure P to its energy density p and captures the thermodynamical properties in a form relevant for 
isotropic space-times in general relativity. Often, one can assume the equation of state parameter w to be constant 
during successive phases of the universe evolution, with sharp jumps between different phases such as w = — 1 during 
^— ^ inflation, followed hy w = during radiation domination and w — 1 during matter domination. 

1 I Observationally relevant details can depend on the precise values of w at a given stage, in particular if one uses 

an effective value describing a mixture of different matter components. For instance, during big bang nucleosynthesis 
one is in a radiation dominated phase mainly described by photons and relativistic fermions. Photons, according to 
Maxwell theory, have an exact equation of state parameter w = i as a consequence of conformal invariance of the 
equations of motion (such that the stress-energy tensor is trace-free). For fermions the general equation of state is 
^ more complicated and non-linear, but can in relativistic regimes be approximately given by the same value w = ^ 
^ [ as for photons. In contrast to the case of Maxwell theory, however, there is no strict symmetry such as conformal 
invariance which would prevent w to take a different value. It is one of the main objectives of the present paper to 
discuss possible corrections to this value. 

For big bang nucleosynthesis, it turns out, the balance between fermions and photons is quite sensitive. In fact, 
different values for the equation of state parameters might even be preferred phenomenologically Jj]. One possible 
reason for different equations of state could be different coupling constants of bosons and fermions to gravity, for 
which currently no underlying mechanism is known. In this paper we will explore the possibility whether quantum 
gravitational corrections to the equations of state can produce sufficiently different values for the equation of state 
parameters. In fact, since the fields are governed by different actions, one generally expects different, though small, 
correction terms which can be of significance in a delicate balance. Note that we are not discussing ordinary quantum 
corrections of quantum fields on a classical background. Those are expected to be similar for fermions and radiation 
in relativistic regimes. We rather deal with quantum gravity corrections in the coupling of the fields to the space-time 
metric, about which much less is known a priori. Thus, different proposals of quantum gravity may differ at this 
stage, providing possible tests. 

An approach where quantum gravitational corrections can be computed is loop quantum cosmology [2], which 
specializes loop quantum gravity 0, Q to cosmological regimes. In such a canonical quantization of gravity. 



2 



equations of state must be computed from matter Hamiltonians rather than covariant stress-energy tensors. Quantum 
corrections to the underlying Hamiltonian then imply corrections in the equation of state. This program was carried 
out for the Maxwell Hamiltonian in and is done here for Dirac fermions. There are several differences between the 
treatment of fermions and other fields, which from the gravitational point of view are mainly related to the fact that 
fermions, in a first order formulation, also couple to torsion and not just the curvature of space-time. After describing 
the classical derivation of equations of state as well as steps of a loop quantization and its correction terms, we use 
big bang nucleosynthesis constraints to see how sensitively we can bound quantum gravity parameters. 

We have aimed to make the paper nearly self-contained and included some of the technical details. Sees. Illll on the 
canonical formulation of fermions, IIVI on quantum corrections from loop quantum gravity and FVl on the analysis of 
big bang nucleosynthesis can, however, be read largely independently of each other by readers only interested in some 
of the aspects covered here. We will start with general remarks on the physics underlying the problem. 

II. THE PHYSICAL SETTING 

Big bang nucleosynthesis happens at energy scales E'bbn ~ MeV which are large, but still tiny compared to the 
Planck energy Mp. Also the universe has already grown large compared to the Planck length £p at this stage, and 
space-time curvature is small. One may thus question why quantum gravity should play any role. There is certainly 
a fine balance required for successful big bang nucleosynthesis, but the expected quantum gravity terms of the order 
E/Mp, obtained based on dimensional arguments, would have no effect. 

However, dimensional arguments do not always work, in particular if more than two parameters Lj of the same 
dimension, or any large dimensionlcss numbers are involved. Then, precise calculations have to be done to determine 
which geometric means Y[i L^' with x/ = may appear as coefRcients, or which powers of dimensionlcss numbers 
occur as factors of correction terms. In loop quantum gravity, we are in such a situation: there is the macroscopic 
length scale L, which in our case we can take as the typical wave length of fields during nucleosynthesis, and also 
the Planck length ip = V Gh which arises due to the presence of Newton's constant G and Planck's constant h. In 
addition, there is a third and in general independent scale £ given by the microscopic size of elementary spatial patches 
in a quantum gravity state. This is a new feature of the fundamentally discrete theory, for which the precise state of 
quantum gravity plays an important role. Although £ must be proportional to the Planck length, its specific value 
for a given state can differ numerically. Then, a detailed calculation must show how L, i and £p appear in quantum 
gravity corrections and which numerical values may arise. 

Alternatively, one can work with only two length scales, L and £p, but one has to deal with a large dimensionlcss 
parameter J\f given by the number of discrete patches of the underlying state in the volume considered, for instance a 
volume of the size such that J\f = L^/£^. Examples of cosmological phenomena are known where this does play a 
role for quantum gravity corrections 0, Q i and here we analyze which features arise in the presence of fermions and 
especially for big bang nucleosynthesis. 

There are precedents where such considerations have played important roles. Best known is the evidence for the 
atomic nature of matter derived by Einstein from the phenomenon of Brownian motion. Also here, there are several 
orders of magnitude between the expected size of molecules and the resolution of microscopes at that time. However, 
there is also a large number of molecules which by their sheer number can and do leave sizeable effects on much larger 
suspended particles. There is, of course, never a guarantee that something analogous has to happen elsewhere. But 
this is to be checked by calculations and cannot always be ruled out based only on dimensional arguments. 

We can use the result of corrections to the equation of state of radiation derived in , but a new analysis is required 
for fermions and their specific action. Even in relativistic regimes, the coupling of fermions to gravity differs from 
other fields, e.g. by torsion contributions which arise already from the kinetic term of the Dirac action. One could 
thus expect that quantum corrections for fermions differ from those to radiation and thus, by throwing off the balance 
during nucleosynthesis, possibly enhance the effect of quantum gravity corrections. Whether or not this happens 
cannot be decided without detailed calculations as they are reported and applied here. 

III. CANONICAL FORMULATION OF DIRAC FERMIONS 

For fermions, one has to use a tetrad rather than a space-time metric y^j,, which are related by e^e^ = g^^, 
in order to formulate an action with the appropriate covariant derivative of fermions. This naturally leads one to 
a first-order formalism of gravity in which the basic configuration variables are a connection 1-form and the tetrad. 
In vacuum the connection would, as a consequence of field equations, be the torsion-free spin connection compatible 
with the tetrad. In the presence of matter fields which couple directly to the connection, such as fermions, this is no 
longer the case and there is torsion |[9|, if^, [T3| ■ 
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A. Non-minimally Coupled Einstein-Dirac Action 



A canonical analysis of gravity minimally coupled to fermions is presented in detail in (l7j . In (l6j , a non-minimally 
coupled action is proposed in order to eliminate violations of parity. Here, we first present a general canonical analysis 
providing a consistent canonical formulation for all coupling parameters. At the end of this section, we briefly comment 
on parity which plays a role in this context. 

The basic configuration variables in a Lagrangian formulation of fermionic field theory are the Dirac bi-spinor 
^ = (■0 ffy'" and its adjoint ^ — (^'*)"^7'^ with 7" being the Minkowski Dirac matrices. The 2-component SL(2,C)- 
spinors ij} and 77 will be used later on in a Hamiltonian decomposition of the action. Being interested in applications 
to highly relativistic regimes, we only deal with massless fermions. Their non-minimal coupling to gravity can then 
be expressed, in the notation of [ll], by the total action 

S'[e,cj,«'] = SG[e,uj] + SY[e,Lo,-^] 
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where a is the parameter for non-minimal coupling and minimal coupling is reproduced for a ^ c». The action 
is composed of the matter contribution S'p resulting from the fermion field and the gravitational contribution Sq 
expressed in terms of 
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where 7 is the Barbero-Immirzi parameter [lOl . [ll| . 

The specific form of the gravitational action is the one given by Hoist , formulated in terms of a tetrad field 
with inverse e^, whose determinant is e. (Thus, the space-time metric is g^i, — e^e^. For all space-time fields, 
/, J, . . . = 0, 1, 2, 3 denote internal Lorentz indices and /i, j/, . . . = 0, 1, 2, 3 space-time indices.) The Lorentz connection 



in this formulation is an additional field independent of the triad, and F^J^{uj) = 2(9[^aj^j"' -I- [u;^,Wy]^'^ is its 



curvature. It also determines the covariant derivative V^^ of Dirac spinors by 

1 



l^"l[ilJ] 



[Vm,V.] = -F^V7[/7j] 



(4) 



in terms of Dirac matrices 7/ (which will always carry an index such that no confusion with the Barbero-Immirzi 
parameter 7 can arise). 



B. Dirac Hamiltonian 



As usually, a Hamiltonian formalism of gravity requires a space-time foliation Stit = const such that one can 
introduce fields and their rates of change, which will provide canonical variables. This is done by referring to a time 
function t as well as a time evolution vector field such that t^^V ^J. ~ 1. Rates of change of spatial fields will then 
be associated with their derivatives along t"' . 

For convenience, one decomposes i'* into normal and tangential parts with respect to by defining the lapse 
function N and the shift vector such that = Nn'-'- + N^^ with N^^n^ = 0. Here, is the unit normal 
vector field to the hypersurfaces Sj. The space-time metric induces a spatial metric qfj,u{i) on St by the formula 
.9^11^ = <l^J^v ~ niinv This is one of the basic fields of a canonical formulation, and its momentum will be related to qab 
defined as the Lie derivative of qab along t°- . Since contractions of q^i, and N^^ with the normal vanish, they give 
rise to spatial tensors qab and N"-. 

In our case, we are using a tetrad formulation, where provides a map from the tangent space of space-time to an 
internal Minkowski space. The space-time foliation thus requires an associated space-time splitting of the Minkowski 
space. This takes the form of a partial gauge fixing on the internal vector fields of the tetrad: the directions (or rather 
boosts) of tetrad fields can no longer be chosen arbitrarily. Instead, we decompose the tetrad into a fixed internal unit 
time-like vector field and a triad on the space Sj. We choose the internal vector field to be constant, n/ = —Sj^ with 
n^n/ = — 1. Then, we allow only those tetrads which are compatible with the fixed in the sense that n° — n^Cj 
is the unit normal to the given foliation. This implies that ej ^ £f — n'^nj with Efua — £fn^ = so that £f is a 
spatial triad. 
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Now, using ej=£f — n°'ni with n/ = — ^/,o and n°' = N "^{t"" — N°-) we can decompose the non-minimally coupled 
Dirac action and write it in terms of spatial fields only: 

(5) 

where the space-time determinant is factorized as e = N ^ with the determinant q of the spatial metric. These terms 
can be decomposed into several terms containing the SL(2, C)-spinors ■0 and 77 instead of the Dirac spinor 
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- (ie'„„a;r"(V'V - V^v) + ie^„„c. ™o J") , (6) 



and similarly an expansion of the terms involving the non-minimal coupling parameter a in the action ^ gives 
la 

+ £ (^V^9a^ + Tl^o'daTl + (aa^/'^)^^ -f (aaTy^)^'?? + C^^'^ Jfc " Jq) , (7) 

where c.c. denotes complex conjugation and we have introduced the fermion current J*^ := ip^a'^tp + rj^cr^rj and the 
time component Jq := f/'^''/' ~ Tj'^r]. Details of this calculation as well as the corresponding canonical decompositions 
of the gravitational part can be found in • 

We also refer to [13] for the definition of canonical gravitational variables as they appear also in the matter action. 
In terms of the connection components w^"^ in Eq. ([¥]), we define 



, Kl:^-u^ (8) 



2 ^ 

and the Ashtekar-Barbero connection [l^, [3 

Al Tl + jKi = ~le\^uj^' - ^u;J° (9) 

with the Barbero-Immirzi parameter 7 [Tlj . This connection is important because it provides a convenient canonical 
structure with Poisson brackets 

{Alix), E'^iy)} = 87T^G5'j}6{x - y) (10) 

with the densitized triad Ej v^^j"- Moreover, following the steps presented in [13], it is straightforward to show 
that the torsion contribution to the spin connection T is given by 

C^^^TTT^ffl-r) ^'u^aJ'-il+^AeiA. (11) 



4(1 + 72) \V a' 

One can then show that = -I- C* , with T\ being the torsion- free spin connection compatible with the co-triad 

For a fixed value of 7, the above equation for torsion reduces to that for minimal coupling as a — > 00. As noted in 
p^ . based on a Lagrangian analysis, another interesting case is a = 7, which reduces to — ■|e^J'^, making torsion 
7-independent. The non-minimal coupling is important for the behavior of the action under parity reversal, and it 
is the value a = 7 which provides parity invariance of the combined system of gravity and fermions. This may seem 
surprising because neither vacuum gravity nor the minimal Dirac action on a fixed background violate parity, and thus 
the extra term of non-minimal coupling seems to introduce parity violation. As [l6j in the Lagrangian picture and 
[17I ] in the Hamiltonian picture show, however, the torsion introduced by coupling fermions to gravity also introduces 
parity violation in the Hoist action, which has to be canceled by the additional term of non-minimal coupling. (The 
action with minimal coupling was called "not fully consistent" in (l6| , while our analysis is consistent for any value of 
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a. To avoid potential confusion, one should first note that flq starts from the Einstein-Cartan action with minimal 
coupling, not directly from a Hoist action as here. Parity properties are different in both cases, and the observation 
of [lal refers only to this behavior. The Hoist action with minimal coupling would indeed be inconsistent with parity 
preservation, but it does not present an inconsistency of the overall framework unless one explicitly requires parity 
preservation. In fact, parity violation is expected at the onset of BBN due to weak interactions. Anyway, the following 
discussion of our paper can remain unchanged for different values of a and thus applies directly to the parity preserving 
case.) 

Upon inserting (fTTj) in ^ and (O, the action (O takes the form 



\ ran jk 



+ \^lknKEl (l - ^) J" 



(12) 



where Oi^ := i(l + ^), 9^ := -1(1 — ^) and we have used the covariant derivatives, 'Da ~ da + A^^ti, related to the 
Ashtekar-Barbero connection, and t; = —^iai. The first part of the action in this form shows that momenta of the 
fermion field are = of ^ = ^^"0 and gv — ~ix^ of x = -^rj. Here, we are using half-densitized spinor fields 
to avoid complex- valued canonical variables [1^ . 

Moreover, the action provides the Hamiltonian H = J (PxN{x)6S/SN. Together with fermion dependent terms 
resulting from the gravitational action, this provides derivative terms and self-interaction terms in the Hamiltonian 



^^Dirac — 



c^xN 



St 



'^Va (ttJt^C + ~ {eLnJr^Va^ - 9nn^r^VaX - c.c.) 



(13) 



7K/3 



with /9 := 7 + ^. The top line of this expression is the most important one because its derivative terms are dominant 
in relativistic regimes. In addition to those, we highlight the presence of four-fermion interactions in the second line, 
which we summarize as 



B := 



(3 - 2 + 27^) {.Jni + .^nxK^fr^C + n^x) + 1^ (i - 2) - -^x)(-[C - ^Ix) (14) 



2(1 + 7^) 
as it multiplies g^^/^ 



C. Equation of state 

From the Hamiltonian we can determine energy and pressure and formulate the equation of state. The matter 
Hamiltonian is directly related to energy density by 



1 ^ilDirac 



and thus, from p3p . the energy density is 



'^^^ ( -^Va {irfr^^ + ttJt^x) + * (-^L^r^^aC + 9R7T^T^VaX - c.c.)] + - (16) 



The canonical formula for pressure is 



P= ?_£;a^^Dirac , . 
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as shown by a straightforward adaptation of the calculation done in fd\ for metric variables. Now using the functional 
derivative 



and thus 



^^ = 2^.^(--y), (18) 



' ^ - —mSi-e'^4)K^-y), (19) 



and inserting (|19p in ^T7\ . we obtain the pressure 



- (20) 



This results in an equation of state 



P 1 2B 

WDirac = " = ^ ^ 1" • (^1) 

p 3 3p 



In relativistic regimes, the kinetic term involving partial derivatives da contained in Va is dominant, which leaves 
us with an equation of state 

P 1 

w = - = - + e (22) 
P 3 

whose leading term agrees with the parameter for a Maxwell field. But there are clearly correction terms for fermions 
already in the classical first order theory. They do not arise for the Maxwell field, implying a difference in the coupling 
to gravity due to torsion, which is present even in relativistic regimes. The order of magnitude of the additional term 
depends on the fermion current density and is thus not expected to be large unless regimes are very dense. We will 
not consider this correction further in this article, but highlight its role as a consequence of torsion. 

In addition, the canonical analysis performed here provides the stage for a quantization of the theory, resulting in 
further correction terms from quantum gravity as they occur even for the Maxwell field 0]. Their magnitudes must 
also be extracted and then compared for the different fields. 



IV. QUANTUM CORRECTIONS 



We are now in a position to derive the form of quantum gravity corrections for fermion fields in an explicit scenario 
of quantum gravity, completing the program of Sec. |TT1 We will not require the quantized dynamics of gravity, which 
makes our aim feasible. But we do need to know how the triad is quantized because it appears in the Dirac Hamiltonian. 
It appears as a basic variable in the Poisson relations (fTO|) . together with the Ashtekar-Barbero connection A^. Loop 
quantum gravity Q, 0] proceeds by forming a quantum representation of these basic fields which turns holonomies 
he{A) = V exp{J^dte''Al^Ti) along arbitrary spatial curves e into creation operators of geometrical excitations along 
these curves. Many such creations then excite a discrete state based on the graph formed by all curves used. On such 
a state, Ef is quantized through the fluxes Fs{E) — d^yEfUaT'^ through surfaces S with co-normal n^. Resulting 
operators have eigenvalues obtained by counting intersections of the surface and the graph in the state, weighed by 
multiplicities of overlapping curves. 

This structure is important for a background independent representation which replaces the usual Fock represen- 
tation of quantum fields on a curved background space-time. The Fock representation requires a metric to define 
basic concepts such as a vacuum state or creation and annihilation operators, which can thus not be used for the 
non-perturbative gravitational field itself. This mathematical property has important physical implications because 
extended objects, namely the integrated holonomies and fluxes, are used as basic quantities. In this way, there is 
non-locality and other implications which affect quantum gravity corrections. 

In loop quantum gravity, there are three main effects which imply correction terms in effective matter equations. Any 
Hamiltonian contains inverse powers of triad components such as in (jl3p . which for fermions is a consequence of 

the fact that they are quantized through half-densities [l^. The loop quantization, however, leads to triad operators 
which have discrete spectra containing zero, thus lacking inverse operators. A proper quantization, along the lines of 
p~9l] . does give a well-defined operator with the correct semiclassical limit. But there are deviations from the classical 
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behavior on small length scales, which are the first source of correction terms. As in the case of the Maxwell field [6|, 
this is the main effect we include here. 

In addition, there are qualitatively different correction terms. First, loop quantum gravity is spatially discrete, 
with states supported on spatial graphs. Quantizations of Hamiltonians thus lead to a discrete representation of any 
spatial derivative term as they also occur for fcrmions. The classical expression arises in a continuum limit, but for 
any given state the discrete representation implies corrections to the classical derivatives as the leading terms in an 
expansion. Secondly, the connection is quantized through holonomies rather than its single components. Thus, the 
quantum Hamiltonians are formulated in terms of exponentials of line integrals of the connection which also give the 
leading classical term plus corrections in an expansion. Finally, whenever a Hamiltonian is not quadratic, there are 
genuine quantum effects as they occur in typical low energy effective actions. They can be computed in a Hamiltonian 
formulation as well [13, HH) contributing yet another source of corrections. 

One certainly needs to know the relative magnitude of all corrections in order to see which ones have to be taken 
into account. For all of them, the magnitude depends on details of the quantum state describing the regime. Here, 
properties of states have to be taken into account, and dimensional arguments are no longer sufficient. For instance, 
discretization and curvature corrections depend on the patch size occurring in the discrete state underlying quantum 
gravity. This patch size is typically small compared to scales on which the matter field changes, even in relativistic 
regimes assumed here. Thus, such corrections can be ignored in a first approximation. What remains are corrections 
from inverse powers. While other corrections shrink in the continuum limit where the patch size becomes small, inverse 
corrections actually grow when the patch size approaches the Planck length. The regimes where the two classes of 
corrections are dominant are thus neatly separated, and we can safely focus on inverse triad corrections only. The 
relevant formulas are collected in Appendix [XJ see also [l^]. A detailed and complete derivation is not yet available 
since precise properties of a quantum gravity state would be required. Still, many general qualitative insights can be 
gained in this way. 

In the Dirac Hamiltonian (|13p the factor to be quantized containing inverse powers of the densitized triad is 

9 pa pj pk e2 3 k 



' ' Vv 

in terms of the volume Vy ~ <l{v) of one discrete patch at a point v. We can already notice the close resemblance 
to the Maxwell Hamiltonian, where the corresponding expression is qab/f-G^/q = (-o^/q which differs only by the 
additional e-tensors. This close relation will, in the end, lead to quite similar quantum corrections for photons and 
fermions. 

We proceed using (|A6p for r = 1/2, and write 

2Ef ( £o 



abc 



e,,feK,K,nMc,V;^}, (23) 



which can then be quantized by turning Poisson brackets into commutators of operators. This results in 

^-^f _ .KIJ^ (P V-^^pJ\(P \r-^pk\ _ pKU ^ ^(1/2)^(1/2) .J .fe 

j^-j^-e eijk{ioVv ej){toVv ej)-e ^ijkty^j C^j c)(/)<3(j) i24J 

with Cl^j defined in (|A7|) . As explicit calculations in the appendix show, eigenvalues (jA9p of the operator C^^^l^\ or 
its expectation values in semiclassical states, do not follow precisely the behavior expected for the classical function it 
quantizes. Deviations become larger for small values of the triad components, which can be captured in a correction 
function multiplying the classical expression. In particular, such a correction function will appear in a 

Hamiltonian operator, thus also correcting expressions for energy density and pressure or the equation of state. The 
functional form follows from Eqs. (|A9|) and (jAlOp . 

Thus the general expression one can expect for a phenomenological Dirac Hamiltonian including corrections from 
inverse powers of the triad is 

r f j?a 9(E'') \ 

ffphcn = / d^xN (~/32?a {ttJt^^ + TT^T^x) - 2z {eLTrJr'Va^ - enn^T^VaX - c.c.)) + ^-^B (25) 

with two possibly different correction functions a and 9. This also affects the energy density and pressure terms, 
derived by the general expressions (|15p and (|17p . We are mainly interested in the correction to the one-third in the 
equation of state (HJ), so we focus on the first term in in what follows. Moreover, for a nearly [i^ isotropic 
background geometry a only depends on the determinant q of the spatial metric and thus q°'^5a/ Sq""^ = —3qda/dq = 
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FIG. 1: The correction function p7|) as a function of the scale factor (sohd line). The asymptotic form (|28|l for large a is shown 
by the dashed line. (The sharp cusp, a consequence of the absolute value appearing in (|27|) . is present only for eigenvalues as 
plotted here, but would disappear for expectation values of the inverse volume operator in coherent states. This cusp will play 
no role in the analysis of this paper.) 



— iflda/da with the scale factor a related to qhy q = det{qab) = ■ In this case the quantum gravitational expectation 
for a{q), as per Eqs. (jA9p and (jAlOp . simplifies. To use these expressions, we have to relate the scale factor to quantum 
gravitational excitation levels as they occur in calculations of loop quantum gravity. In the notation of the appendix, 
an elementary discrete patch in a nearly isotropic space-time has, on the one hand, an area of ^g*^^ if is the 
coordinate diameter of the patch. This can be expressed as ^pa^ = (Vy/A/y)^^"^ where A/y is the number of patches 
in a box V of volume Vy. On the other hand, using (jA2[) the quantum gravity state assigns a value of An'-fip^y to 
this patch via the flux operator, where fXy is the quantum number of the geometrical excitation of this patch. Thus, 
we obtain 



where 



1/3 



adisc = 2^ip ] (26) 



Vo 



with the coordinate volume Vq of the box V. The numerical value of adisc depends on coordinates via Vq, or on the 
normalization of the scale factor. (It does not depend on the choice of the box V because a change would multiply 
A/y and Vo by the same factor.) But it is important to note that adisc is not just determined by the Planck length £p, 
which appears for dimensional reasons, but also depends on the large number A/y of discrete patches per volume as 
given by the quantum gravity state. This is exactly a parameter as expected in the discussion of Sec. |TT1 Replacing 
fii, in the equations of the appendix, we obtain 

a(a) = 8V2(a/adisc)' ((2(a/adisc)' + 1)'/' - |2(a/adisc)' - 11'/')' (27) 

where adisc appears, influencing the size of quantum gravity corrections. 

The function is plotted in Fig.[T] One can easily see that a{a) approaches the classical value a = 1 for a 3> adisc/ V2, 
while it differs from one for small a. For a > adisc/\/2, the corrections are perturbative in a^^. 
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This is the first correction in an asymptotic expansion for eigenvalues. If semiclassical states rather than volume 
eigenstates are used, powers of in the leading corrections can be smaller. Moreover, via Afy the discreteness scale 
adisc is expected to be not precisely constant but a function of a itself because the underlying spatial discreteness of 
quantum gravity can be refined dynamically during cosmological evolution l^HSl- (Indeed, dynamical refinement is 
also required for several other phenomenological reasons [25,, M M M In 

our following analysis we will thus 

assume a functional form 

a(a) = 1 + c(a/ao)"" (29) 

where we traded the fundamental normalization by Odisc for normalization with respect to the present-day value of 
the scale factor ag. From the derivation, n is likely to be a small, even integer and c is known to be positive. The 
constant c depends on adisc and inherits the A/y-factor. It can thus be larger than of order one. We will treat this 
parameter as phenomenological and in the end formulate bounds on c as bounds for A/y- 
Energy density and the pressure then are, ignoring the classical interaction term B, 

Pes = ^a(a) {-(3Va {nfr'^ + ttJt^x) + i (-^i^f r'PaC - 9Rn^T'VaX - c.c.)) (30) 



and 



3Peff = ^a(a) (l - i-PV, {irlr^i + tt^t'x) + * (-^LTrJr'P.e - 0«<r'I?.X - c.c.)) . (31) 



From this, the equation of state w can easily be computed: 



1 / dloga' 



3 \ d log a 

This quantum gravity correction is independent of the specific matter dynamics as in the classical relativistic case. 
It results in an equation of state which is linear in p, but depends on the geometrical scales (and the Planck length) 
through a. This is the same general formula derived in Q for radiation. Thus, on an isotropic background radiation 
and relativistic fermions are not distinguished by the form of quantum corrections they receive. 
With the equation of state and the continuity equation 

/3 + 3-(p + P) = 0, (33) 
a 

where a is the scale factor and the dot indicates a proper time derivative, we can determine energy density as a 
function of the universe size. We first obtain 

dlogp(a) 

— = -3 1 + wa) , (34 

dloga 

explicitly showing the dependence of the equation of state on the scale factor. The solution is 



p{a) = po exp 



-3 / (l + w(a)) dloga 



(35) 



where po is the integration constant. Inserting an equation of state of the form (15^ we obtain 

p{a) ~ pQa{a)a^^ . (36) 

For a = 1, we retrieve the classical result p{a) oc a~^, but for a ^ 1 loop quantum gravity corrections induced by 
the discreteness of flux operators are reflected in the evolution of energy density in a Friedmann-Robertson- Walker 
universe. 

V. EFFECT ON BIG BANG NUCLEOSYNTHESIS 

The production of elements in the early universe is highly sensitive to the expansion rate, given by 

1/2 



, x^Gp , (37) 
a V 3 ' 



10 



where p is the total density, thus including radiation and fermions. As we have seen here for fermions and in Q 
for radiation, the effect of loop quantum gravity corrections is to multiply the effective p{a) by a factor a{a). Most 
importantly, we find that a (a) is the same for both bosons and fermions (up to possible quantization ambiguities), so 
a separate treatment of the two types of particles in the early universe (as in Ref. [H) is unnecessary here. 

In the standard treatment of the thermal history of the universe, the density of relativistic particles (bosons or 
fermions) is given by 

P = ^9*T\ (38) 

where 5* is the number of spin degrees of freedom for bosons, and 7/8 times the number of spin degrees of freedom 
for fermions, and T is the temperature, which scales as 

Toca^^ (39) 

The equation of state parameter is 

w = 1/3. (40) 

Clearly, equations ([38|) - ((40|) are inconsistent with equations (|32p and ((36|) . There is some ambiguity in determining 
the correct way to modify the expressions for p{T) and w. We have chosen to assume that the modifications are 
contained in the gravitational sector, so that the density is given by 

P = o^{af^9*T\ (41) 



with the temperature scaling as in equation (|39|) . and the equation of state w is given by equation (j32|l . This 
guarantees that the standard continuity equation (j33p continues to hold. Note that this is not the only way to 
incorporate equation (|36|1 into the calculation, but it seems to us the most reasonable way. This issue requires a 
consideration of thermodynamics on a quantum space-time, which is a fascinating but not well-studied area. Instead 
of entering details here, we note that we interpret the a-corrcction as a consequence of a quantum gravity sink to 
energy and entropy. Thus, quantum gravity implies a non-equilibrium situation which would otherwise imply that p 
must be proportional to T'* without any additional dependence on a oc T^^ . 

With these assumptions, we can simply treat a{a) as an effective multiplicative change in the overall value of G. 
Note that this simplification is only possible because we explicitly derived by our canonical analysis that, unexpectedly, 
quantum corrections of radiation and fermions appear in similar forms. This makes possible a comprehensive derivation 
of implications for BBN, bearing on earlier work. In fact, a great deal of work has been done on the use of BBN to 
constrain changes in G (see, e.g. Refs. [13, HH, HI, [H, [sj] ) . The calculation is straightforward, if one has a functional 
form for the time-variation in G. For the loop quantum gravity corrections considered here, the most reasonable 
functional form is (|29p . Note that this expression is by construction valid only in the limit where a{a) — 1 << 1. In 
terms of the effective gravitational constant, G, one can then write 

G(a) = Go[l + c(a/ao)""], (42) 

where Go is the present-day value of the gravitational constant. 

In order to constrain the values of c and n, we calculate the predicted element abundances with the indicated 
change in G and compare with observational constraints. Big Bang nucleosynthesis proceeds first through the weak 
interactions that interconvert protons and neutrons: 

n + Ve ^ p + , 

n ^ p + e^ + Ve- (43) 

When T ^ 1 MeV, the weak-interaction rates are faster than the expansion rate, a/a, and the neutron-to-proton ratio 
(n/p) tracks its equilibrium value exp[— Am/T], where Am is the neutron-proton mass difference. As the universe 
expands and cools, the expansion rate becomes too fast for weak interactions to maintain weak equilibrium and n/p 
freezes out. Nearly all the neutrons which survive this freeze-out are converted into '^He as soon as deuterium becomes 
stable against photodisintegration, but trace amounts of other elements are produced, particularly deuterium and ^Li 
(see, e.g., Ref. [35| for a review). 

In the standard model, the predicted abundances of all of these light elements are a function of the baryon-photon 
ratio, rj, but any change in G alters these predictions. Prior to the era of precision CMB observations (i.e., before 
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WMAP), Big Bang nucleosynthesis provided the most stringent constraints on ij, and modifications to the standard 
model could be ruled out only if no value of rj gave predictions for the light element abundances consistent with the 
observations. However, the CMB observations now provide an independent estimate for 77, which can be used as an 
input parameter for Big-Bang nucleosynthesis calculations. 

Copi et al. 33] have recently argued that the most reliable constraints on changes in G can be derived by using 
the WMAP values for 77 in conjuction with deuterium observations. The reason is that deuterium can be observed in 
(presumably unprocessed) high-redshift quasi stellar object (QSO) absorption line systems (see Ref. ^S^l and references 
therein), while the estimated primordial *He abundance, derived from observations of low metallicity HII regions, is 
more uncertain (see, for example, the discussion in Ref. |13|)- While we agree with the argument of Copi et al. in 
principle, for the particular model under consideration here it makes more sense to use limits on '^He than on deuterium, 
in conjunction with the WMAP value for ij. The reason is that the ^He abundance is most sensitive to changes in 
the expansion rate at T ^ 1 MeV, when the freeze-out of the weak interactions determines the fraction of neutrons 
that will eventually be incorporated into ^He. Deuterium, in contrast, is produced in Big Bang nucleosynthesis only 
because the expansion of the universe prevents all of the deuterium from being fused into heavier elements. Thus, 
the deuterium abundance is most sensitive to the expansion rate at the epoch when this fusion process operates 
(T ~ 0.1 MeV). The importance of this distinction with regard to modifications of the standard model was first noted 
in Ref. [3^, and a very nice quantitative analysis was given recently in Ref. [1^. Note that our estimate for the 
behavior of G{a)/Go — 1, equation (|42p . is a steeply decreasing function of a. Thus, the change in the primordial ''He 
abundance will always be much larger than the change in the deuterium abundance. Therefore, we can obtain better 
constraints on this model by using extremely conservative limits on ^He, rather than by using the more reliable limits 
on the deuterium abundance. For the same reason, we can ignore any effect on the CMB, since the latter is generated 
at a much larger value of a, and any change will be minuscule. Hence, we can confidently use the WMAP value for ry. 

WMAP gives 

V = a.iwt^ofZ X 10^'". (44) 

Because the estimated errors on 77 are so small, we simply use the central value for 77; the bounds we derive on c in 
equation (|29p change only slightly when 77 is varied within the range given by equation (j44p . Since c in equations (j29p 
and (|42l) is thought to be positive, the effect of LQG corrections is to increase the primordial expansion rate, which 
increases the predicted ^He abundance. We therefore require an observational upper bound on the primordial "^He 
abundance. As noted earlier, this is a matter of some controversy. We therefore adopt the very conservative upper 
bound recommended by Olive and Skillman [37[ : 

Yp < 0.258, (45) 

where Yp is the primordial mass fraction of *He. For a fixed value of n in equation (H^ . we determine the largest 
value of c that yields a primordial ^He abundance consistent with this upper limit on Yp. Since we are essentially 
bounding the change in G at a/ao ~ 10^^", it is convenient to rewrite equation ([29]) as 

a^l + c/a'^o (46) 

where aio = 10^*'(a/ao). This upper bound on c as a function of 77, is given in Fig. [2] For the special case 77 = 4, we 



can use these results to place a bound on Odisc in equation (|28p . We obtain 

^< 2.4x10-1°. (47) 
ao 

This is not a strong bound for the parameters of quantum gravity, but clearly demonstrates that quantum corrections 
are consistent with successful big bang nucleosynthesis. 

In terms of more tangible quantum gravity parameters, we have 

K" < m 



for the number of patches at the time of big bang nucleosynthesis. In terms of the volume V\ = (10-i°ao)'^Vb at this 
time, we have A/y < 3Vv/^p with the value 7 k, 0.24 of the Barbero-Immirzi parameter as derived from black hole 
entropy calculations. More meaningfully, if we view ^^Fiv^tp as the basic length scale as it appears in the spectrum 



()A2[) of loop quantum gravity, the bound becomes more interesting: This gives A/'y^'^ < 2.W^^'^ ji^^fiFilp). This 



upper limit is already quite close to what one expects for elementary patch sizes in loop quantum gravity, which 
would provide A^y^ < Y^''^ l{2^ip) as a fundamental upper limit. Given that these values are close to each other, 
we see a clear potential of improvements by more precise observational inputs. Moreover, other correction terms from 
quantum gravity could be used to obtain a lower bound for Ay such that the allowed window would be reduced to a 
smaller size. 
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FIG. 2: Solid curve gives upper bound on c as a function of n, for the assumed form for a: a — 1 + c/aiQ, where aio is the 
value of the scale factor in units for which aio = 10^" at present. 

VI. CONCLUSIONS 

Big bang nucleosynthesis is a highly relativistic regime which, to a good approximation, implies identical equations 
of state for fermions and photons. There are, however, corrections to the simple equation of state w = ^ for fermions 
even classically. One observation made here is that the interaction term derived in leads to such a correction and 
might be more constrained by nucleosynthesis than through standard particle experiments [Tsj . We have not analyzed 
this further here because more details of the behavior of the fermion current would be required. 

A second source of corrections arises from quantum gravity. Remarkably, while quantum gravity effects on an 
isotropic background do correct the equations of state, they do so equally for photons and relativistic fermions. 
Initially, this is not expected for both types of fields due to their very different actions. Thus, quantum gravity effects 
do not spoil the detailed balance required for the scenario to work and bounds from big bang nucleosynthesis obtained 
so far are not strong. But there are interesting limits for the primary parameter, the patch size of a quantum gravity 
state. It is dimensionally expected to be proportional to the Planck length £p but could be larger. In fact, current 
bounds derived here already rule out a patch size of exactly the elementary allowed value in loop quantum gravity. 
With more precise estimates, these bounds can be improved further. 

We have made use of quantum gravity corrections in a form which does not distinguish fermions from radiation. 
Although the most natural implementation, quite unexpectedly, provides equal corrections as shown here, there are 
several possibilities for differences which suggest several further investigations. Small deviations in the equations of 
state and thus energy densities of fermions and radiation are possible. First, there are always quantization ambiguities, 
and so far we tacitly assumed that the same basic quantization choice is made for the Maxwell and Dirac Hamiltonians. 
Such ambiguity parameters can be exphcitly included in specific formulas for correction functions; see e.g. (2^.l4ll[4^. 
Independent consistency conditions for the quantization may at some point require one to use different quantizations 
for both types of fields, resulting in different quantum corrections and different energy densities. Such conditions 
can be derived from an analysis of anomaly-freedom of the Maxwell field and fermions coupled to gravity, which is 
currently in progress. As shown here, if this is the case it will become testable in scenarios sensitive to the behavior 
of energy density such as big bang nucleosynthesis. Moreover, assuming the same quantization parameters leads to 
identical quantum corrections for photons and fermions only on isotropic backgrounds. Small-scale anisotropics have 
different effects on both types of fields and can thus also be probed through their implications on the equation of 
state. 

For this, it will be important to estimate more precisely the typical size of corrections, which is not easy since it 
requires details of the quantum state of geometry. The crucial ingredient is again the patch size of underlying lattice 
states. On the other hand, taking a phenomenological point of view allows one to estimate ranges for patch sizes 
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which would leave one in agreement with big bang nucleosynthesis constraints. Interestingly, corrections studied here 
provide upper bounds to the patch size, and other corrections from quantum gravity are expected to result in lower 
bounds. A finite window thus results, which can be shrunk with future improvements in observations. 
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APPENDIX A: ELEMENTS OF LOOP QUANTUM GRAVITY 

We collect the basic formulas required to compute quantum corrections in loop quantum cosmology, referring to 
, [IS] for further details. We can restrict our attention to perturbative regimes around a spatially flat isotropic 
solution where one can choose the canonical variables to be given by functions {p^ (x), kj{x)) which determine a 
densitized triad by Ef = p'-^\x)Sf and extrinsic curvature by Kl = This diagonalization implies strong 

simplifications in the explicit evaluation of formulas [H, . 

Any state is associated with a spatial lattice with U(l) elements rjyj attached to its links Cyj starting at a 
vertex v and pointing in a direction along a translation generator Xf of the homogeneous background. In our 
context of perturbations around Friedmann-Robertson-Walker space-times, this can be seen as a definition of the 
quantum gravity state on which the quantized and act. The U(l) elements rjy j appear as matrix elements in 
SU(2) holonomies h^ j — Herjy j + 2Tjliiir]y j and represent the connection, while fluxes F^ j are used as independent 
variables for the momenta. An orthonormal basis of the quantum Hilbert space is given by functionals | . . . , /iu,/, ■ • •) = 
Yiv iVv"/ ^ possible choices of fiyj E Z. Together with basic operators, which are represented as holonomies 

Vvjl ■■■,Hv',j,---) = \---, ^J.vj + 1, • ■ •) (Al) 

for each pair (w, /), and fluxes 

Tvj\ . . . , ^J.v',J, ■ ■ ■) = 2T:-/£p{fj,yj + fJ-v-i)\ ■ ■ ■ ,tJ-v',j, . . .) , (A2) 

this defines the basic quantum representation. Here, ip — \fhG is the Planck length and a subscript — / means that 
the edge preceding the vertex v in the chosen orientation is taken. 

For quantum corrections to classical equations, one important step in addition to computing expectation values of 
operators is the introduction of a continuum limit. This relates holonomies 

j^^j = exp(i / At-f~ki/2) w exp(i47fc/(u + //2)/2) (A3) 

to continuum fields ki through mid-point evaluation (denoted by w -f- 1/2), and similarly for fluxes 

F.,i = I P\v)A^y « IIp\v + 1/2) . (A4) 

Here, Ho is the coordinate length of lattice links, which enters the continuum approximation since we are integrating 

classical fields in holonomies and fluxes. From flux operators one defines the volume operator V — n/=i \J \^v,i\, 

using the classical expression V — J d'^Xy/\p^p'^p'^\ « J2v ^o\/W^W\ — J2v \/\P^P^P^\- With (|A2p . its eigenvalues 
are 

viifi^j}) = {2njelf'Y.f[ + (A5) 

V 7=1 

The volume operator is central for inverse triad corrections because inverse densitized triads, or a co-triad, can be 
quantized using relations such as 

{A^,y;} = 4^7Gry:-^el (A6) 
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for < r < 2. Even though an inverse power of volume, together with a co-trad, occurs on the right hand side, the 
left hand side can be quantized directly in terms of the volume operator, using holonomies for connection components, 
and turning the Poisson bracket into a commutator [l^, [4^1 ■ Such a quantization leads to 

Vf^e^: = E ^^<r^K.AK!..yi\) = " B^l^^,, IcJ].^,, (A7) 

(Te{±i} 

where, for symmetry, we use both edges touching the vertex v along direction X*^ and j is, after taking the trace 
in (IX7| . 



B^vl ■■= 7— TT^Z (sv,iKc,,i - c^jKs^j] (A8) 

with 



Cv,i = \{riv,i + vl,!) and s„,7 = :^(77t-,/ - vlj) ■ 
I ' 2,1 ' 

As in [l^l effects of the quantization of triad (metric) coefficients are included by inserting correction functions in 
the classical Hamiltonian which follow, e.g., from the eigenvalues [13] 

(A9) 

(where indices J and K are defined such that £/ ja- 7^ 0) of operators C^Y^'' ■ 

Classically, we expect qij / ^ = \/\p^P^P^\/p^P'^ for this quantity, with a densitized triad Ef = p^^^Sf and using 
the relation (|A2[) between labels and flux components. Although for large fiyj the eigenvalues indeed approach the 
function 



ciT\{^^.'J'})cif\{^^.'J'}) - (2^74) 



•2^-1/2. 



11^=1 VlfJ'V.K + fJ.v.-K\ 



they from the classical expectation differ for values of iiyj closer to one. This deviation can, for an isotropic back- 
ground, be captured in a single correction function 



av,K-i:2^C^j [{^lv'J'\) --=3 , ==- (AlU) 

which would equal one in the absence of quantum corrections. This is indeed approached in the limit where all 
fJ-vj ^ 1, but for any finite values there are corrections. If all fiyj > 1 one can directly check that corrections are 
positive, i.e. a^.A > 1 in this regime. Expressing the labels in terms of the densitized triad through fluxes (|A2[) results 
in functionals 

a[p\v)] = a^,K(47r74^„j) (All) 

which enter quantum corrections. 
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